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ABSTRACT 

We present the first X-ray detection of resonant absorption from warm/hot 
local gas either in our Galaxy, or in the Intergalactic space surrounding our 
Galaxy, along the line of sight toward the blazar PKS 2155-304. The Chandra 
HRCS-LETG spectrum of this z=0.116 source clearly shows, at > 5cr level, un- 
resolved (FWHM < 800 km s'^ at a 2cr confidence level) OVIIi^^. and NelX^-a 
resonant absorption lines at 21.6031°;°^^ A and 13.448l°:°^2 A (i.e. cz = (UtHo) 
km s~^ in the rest frame, from the OVllxa line). OVlllKa and OVIIi^-^ from 
the same system are also detected at a lower significance level (i.e. ~ 3a), while 
upper limits are set on OVIUkis, NeXxa and NelX/^/j. The FUSE spectrum of 
this source shows complex OVl2s^2p absorption at the same redshift as the X- 
ray system, made by at least two components: one relatively narrow (FWHM 
= 106 ± 9 km s~^) and slightly redshifted {cz = 36 ± 6 km s~^), and one broader 
(FWHM = 158 ± 26 km s'^) and blueshifted {cz = -135 ± 14 km s"^). We 
demonstrate that the physical states of the UV and X-ray absorbers are hard 
to reconcile with a single, purely coUisionally ionized, equilibrium plasma. We 
propose, instead, that the X-ray and, at least the broader and blueshifted UV 
absorber are produced in a low density intergalactic plasma, collapsing towards 
our Galaxy, consistent with the predictions of a Warm-Hot Intergalactic Medium 
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(WHIM) from numerical simulations. We find that any reasonable solution re- 
quires overabundance of Ne compared to O by a factor of ~ 2, with respect to 
the solar value. We propose several scenarios to account for this observation. 



1. Introduction 

The majority of the total baryonic matter in the local {zt^) Universe is predicted to 
be concentrated in highly ionized gas. Structures that are already virialized contain warm 
(10^ — 10^ K) or hot (10^ K) gas [the dense interstellar medium (ISM) of galaxies, and 
the intracluster medium (ICM) of clusters of galaxies]. The greater amount of baryonic 
matter is predicted to lie in, as yet unvirialized, matter in the form of a tenuous warm- 
hot intergalactic medium (WHIM, Hellsten et al., 1998). The detection and study of these 
components is needed for the proper understanding of large and small scale structures in 
the Universe. However, while X-ray emission from the virialized density peaks of the ICM 
and ISM has been detected and intensively studied in X-rays, the predicted highly ionized 
gas in the WHIM has been poorly studied so far, due to instrumental hmitations. The low 
density of the WHIM leads to low emissivity, so that studies of the WHIM in emission are 
a formidable challenge. However, absorption depends only on the total column density of 
the medium, not on density, and background light sources in the form of quasars (Aldcroft 
et al., 1994) and gamma-ray bursts (Fiore et al. 2000) are readily available. A few high 
ionization transitions, notably OVI A1031.93, lie in the UV, but the most prominent ions 
(CVI, OVII, OVIII, NelX) have their strongest transitions in the soft X-ray band (10-40 A) 
for a wide range of temperatures (10^ — 10^'^ K, Nicastro et al., 1999), which should give 
rise to an "X-ray Forest" of absorption hues (e.g. Perna & Loeb, 1998, Fang & Canizares, 
2000). The advent of high resolution {R ~1000) soft X-ray spectroscopy with Chandra and 
XMM-A'"eu;ton allows sensitive studies of the WHIM possible. 

Interstellar OVI was first detected in the UV with the Copernicus sateUite (Jenkins, 
1978a, 1978b, and York, 1977), but only recently has data from the Far Ultraviolet Spec- 
troscopic Explorer (FUSE) sateUite shown the presence of extragalactic OVI intervening 
absorpion (e.g. Tripp et al., 2001, Sembach et al., 2000). 

In this paper we present the first X-ray detection of highly ionized absorption along 
the line of sight towards the blazar PKS 2155-304, with the Chandra HRCS-LETG {High 
Resolution Camera Spectrometer- Low Energy Transmission Grating: Brinkman et al., 2000). 
We argue that the physical and dynamical conditions implied by these lines require an 
extragalactic, low density, origin. Hence these lines are the first detection of the X-ray 
Forest. 
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2. Spectral Fitting 

The bright z=0.116 blazar PKS 2155-304 (Falomo, Pesce, & Treves, 1993) is a caUbration 
source for both the high energy and low energy resohition gratings of Chandra, and has been 
observed several times with all the possible grating-detector configurations. The longest of 
these observation lasted ~ 64.7 ks, for a total net good-time exposure of 62.7 ks, and was 
performed with the HRCS-LETG configuration on 1999 December 25. In the public Chandra 
data archive ^ we found two more HRCS-LETG observations of this source, performed on 
2000 May 31, and 2001 April 6, with net exposures of 25.8 ks and 26.6 ks respectively. 

We retrieved the primary and secondary data products (Fabbiano et al., 2001, in prep.) 
of these three Chandra HRCS-LETG observations of PKS 2155-304, and reprocessed their 
event files with the CIAO (Elvis et al., 2001, in prep.) software (v. 2.1.3) using the most up- 
to-date calibration files as of 2001 August (CALDB v. 6.2) to extract source and background 
spectra and corresponding Ist-order Ancillary Response Files (ARFs) and Redistribution 
Matrices (RMFs) according to the on-line data analysis "threads" provided by the Chandra 
X-ray Center (CXC) ^ . 

During the long 1999 December observation the source was at the particularly high flux 
level of 21 mCrab, with a 0.3-6 keV flux of F ■> = 4.2 x 10"^° erg s~^ cm"^. In both 
the subsequent HRCS-LETG observations, instead, the flux level of PKS 2155-304 was more 
normal at about 6 mCrab in the 0.3-6 keV band. However, despite the large change in 
flux, the spectral shape of the underlying 0.3-6 keV continuum did not change significantly 
between or during the three observations (using simple power law continuum models, we 
measured variations in the spectral slope to be less than AF = 0.1 at 90 % confidence level). 
To increase the signal to noise ratio of our data, we then coadded the three HRCS-LETG 
spectra and averaged the corresponding ARFs, weighting them by the respective exposure 
times. The result (Figure la) is one of the highest quality Chandra grating spectra so far 
for an extragalactic source, with some 600 counts, or a signal to noise ratio of about 25, 
per resolution element (R~ 450, FWHM ~ 660 km s~^ at 20 A) in the continuum, over the 
entire 10-25 A wavelength (~ 0.5 — 1.2 keV) range. 

To ensure a self-consistent analysis, and because of the ~ 30 times higher spectral 
resolution provided, we also retrieved the FUSE (Moos et al, 2000) data of PKS 2155-304 
from the FUSE public archive, extracted the fluxed spectrum, and analyzed the wavelength 
range around OVI (1027-1035 A) and 01 (1037-1040 A) to obtain independent estimates 



^http://asc. harvard.edu/cda/ 

^http://asc. harvard.edu/ciao/docunients_threads. html 
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of the positions, widths and equivalent widths of the local OVl2s->2p(A = 1031.9261) and 
Ol35^3p(A = 1039.23) components, recently discovered and published by Savage et al. (2000) 
and Sembach et al. (2000). The FUSE observation was taken with the LWRS aperture 
(30" X 30") with a total exposure time of 38.6 ks, and covers the entire 980-1180 A range. 
We used data from the LiF mirror and the Al detector segment, which have optimal efficiency 
in the waveband around the OVI line. The data were calibrated and cleaned following the 
procedures presented in the FUSE Data Analysis Cookbook (v. 1.0, Sankrit et al., 2001). 
The resolution of the final data is ~ 20 km s~^. 

We performed simultaneous spectral fitting of both the fluxed FUSE data and the 
responses-folded Chandra HRCS-LETG data, using the Sherpa (Siemiginowska et al., 2001, 
in prep.) package for spectral fitting in CI AO. In the following two sections we summarize 
our results. 



2.1. X-ray Spectral Fitting 

For sources of a few mCrab the inherent order confusion problem of the Chandra HRCS- 
LETG (Brinkman et al., 2000) is important only at A;;25 A. Furthermore, narrow features, 
such as absorption lines, are not in any case seriously affected. So we limited our spectral 
analysis to the shorter wavelengths of 11-24 A (~ 0.5 — 1.2 keV) range, allowing us to use 
only the Ist-order RMFs and ARFs, to fit our spectra. Errors throughout the paper are 
quoted at a 90 % confidence level, for 1 interesting parameter. 

To model the 11-24 A portions of the continuum of PKS 2155-304, we used a power law 
attenuated by an equivalent hydrogen column density of Nh = 1.65 x 10^° cm~^ neutral gas 
to account for the measured Galactic absorption from cold gas along the fine of sight to this 
source (Dickey & Lockman, 1990) ^ . The model used for the neutral absorber includes only 
photoelectric edges by neutral elements, and so does not account for resonant absorption 
lines from the same elements. We obtained a statistically acceptable fit with a power law 
photon index of F = 2.42 ± 0.02, consistent previous results (e.g. Kataoka et al., 2000). 

The residuals to this fit clearly show (at a significance level greater than 8 a) two strong 
and unresolved (see §3.1) features in absorption (Fig. lb,d), at wavelengths consistent with 
the OVIIi^-Q s-nd the NelXi^-Q, resonant lines at zero velocity, plus a third intense 0Iis_»2p 



■^Due to the high clcgrco of saturation of the Galactic disk Nh, the local damped Lya absorption line in 
the HST-STIS spectrum of PKS 2155-304 cannot provide a better estimate of the Galactic column of neutral 
H along this line of sight. 
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resonant absorption line from gaseous 01 (Fig. Ic). Negative residuals are also present 
at the energy of the instrumental resonant molecular 01 absorption line (Fig. Ic). The 
residuals around the expected OVlllKa and the OVII^/? rest-frame wavelengths also show 
absorption features at a lower (~ 3a) significance level (Fig. Ic, see §3.1). We added six 
negative gaussians to our continuum model, and refitted the data, leaving the positions and 
the total flux of all the gaussians free to vary independently, but linking all the widths (in 
units of km s~^) of the high ionization resonance lines from O and Ne. 

The best fitting line parameters values are listed in Table 1, along with the derived ion 
column densities for all the atomic transitions. We also include in table 1 2a upper limits on 
the equivalent widths of the two lines of NeX;^^ and NelX;^^, which we derived by adding 
two more gaussians to our model, freezing their positions to the rest-frame wavelengths of 
these transitions and calculating the errors on the lines' normalizations. 

The strong solid-state instrumental 01 line at ~ 23.3 A (Fig. la,c), is only partly 
accounted for by the current effective area model (Fig. Ic, dashed curve; see also Table 1) 
and is easily separated by the HRCS-LETG from the gas phase 0Iis_>2p resonant transition 
at X{OI) = 23.489 A (K rause, 2001), which is also clearly visible in our data (Fig. Ic), and 
will be object of study of a forthcoming paper. Here we stress only how the intervening 
Galactic cold absorber is clearly distinguished from the highly ionized intervening medium 
seen through the He- like and H-like resonant lines from O and Ne (Table 1; see also §3.1 for 
differences in the dynamics). 



2.2. UV Spectral Fitting 

We also used Sherpa to fit the 1027-1035 A portion of the FUSE spectrum of PKS 2155- 
304 around the strongest line of the OVl2s^2p doublet, at A = 1031.9261 A. Residuals after 
a fit with a power law continuum model, show resolved OVl2s^2p absorption in these data, 
with at least two components present (Figure 2). We added two negative gaussians to our 
best fit continuum model. The fit shows that one line is relatively narrow (FWHM = 106 ±9 
km s^^) and slightly redshifted {cz = 36 ± 6 km s^^ in the rest frame), and the other is 
broader (FWHM = 158 ±26 km s-^) and blueshifted {cz = -135 ± 14 km s'^ see Table 1). 
These values are consistent with those reported by Savage et al. (2000) and Sembach et al. 
(2000). 

We looked for Lya absorption associated with these two components, but 
verified that they both would fall well within the Voigt wings of the damped Lya 
due to cold absorption in our own Galctxy, and eire so not visible. 
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Fig. 1. — Top: (a) X-ray 11-24 A unbinned (AA = 12.5 mA) HRCS-LETG spectrum of PKS 2155-304. 
The folded best fitting model (power-law plus Galactic absorption, plus 6 absorption lines) is shown as a 
solid line. The strongest absorption lines arc identified. Note that the effective area of the HRCS-LETG has 
not been corrected for in this plot, in order to show the strong, but large scale, instrumental features. The 
fitted power-law follows these peaks and dips closely. Bottom: blow-ups of four portions of the spectrum 
in Fig. la: (b) OVIIkq (top left); (c) atomic and sohd-state (i.e. molecular) 01 (top right); (d) NelXi^a 
(bottom left); (c) OWIkh and OVIIIxc (bottom right). Vertical lines indicate the rest frame wavelengths 
of the transitions. The solid-state 01 line in Fig. Ic is an effective area feature (dashed line). 
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Table 1: Best fitting UV and X-ray absorption line parameters 





Line Id. 


A 

(A) 


cz 

(km s"^-) 


FWHM 

(km s--^) 


EW(abs)" 

(cV) 


(10^-" (iir^) 


FUSE 


fWfjN arrow 

ovif™«4 


1032.05 ±0.02 
1031.46 ±0.05 


36 ±6 
-135 ±14 


106 ±9 
158 ± 26 


(2.1 ±0.2) X 10-3 
(1.6 ±0.4) X 10-3 


0.14 ±0.10 
0.11 ±0.30 


Chandra 


OVIIk„ 

OVIIIxa 

OVIIk/3 
NeIXx„ 

NeXxa 
NeIX^f;3 


21.603tn^l 
18.973t°;°^^ 
18.605t°;°i 
13.448t°;°^l 
12.134^* 
11.5467'' 


14+194 

^n+348 

90+490 
^"^-536 

0" 
0" 


< 800 
= FWHM(OVllKa) 
= FWHM(OVIIxa) 
= FWHM(OVIIif„) 
= FWHM(OVIlKa) 
= FWHM(OVllKa) 


(0.311^:?°) 

(0.241°:?°) 
(0.22t°:^6) 
(0.53t°:l) 

< 0.32 

< 0.24 


4.0t2? 

r 9+4.3 

131}^ 

< 5.2 

< 20 




OI(A = 23.489'=) 
Ol(molecular) 


23.509t^;™« 

oq o/i 1+0.026 


-140+114 
NA 


< 170 

< 505 


-(0.251^;}^) 

-(0.301°:?^) 


21 ±10 



" Absorption line EW are formally negative; here we drop the sign. '' Frozen. Krause, 1994. 
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Fig. 2. — Ultraviolet 1027-1035 A fluxed FUSE spectrum of PKS 2155-304, and best fitting model. 
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The FUSE spectrum of PKS 2155-304 also shows the presence of strong Ol3s_»3p(A = 
1039.23) absorption (see Figure 1 in Sembach et al., 2000). For this line we measure A = 
1039.19 ± 0.02 {cz = -11 ± 6 km s"^) and FWHM of 47 ± 5 km s-^ much narrower than 
the width of either of the OVIss^sp components. 

3. Analysis of Absorption Lines 

3.1. Gas Dynamics 

To intercompare the dynamics of the gas producing the detected absorption lines, figure 
3 shows their residuals, in velocity space, from the FUSE and Chandra spectra to their 
best fitting continuum models. Zero velocities in this diagram correspond to the rest frame 
wavelengths of the detected transitions (Verner et al., 1996, Krause, 1994). The spectral 
resolution of the HRCS-LETG (R~ 450, FWHM ~ 660 km s'^ at 20 A) does not resolve 
the high ionization X-ray absorption lines (see Table 1). 

The spectral resolution of the FUSE spectrometer {R ~ 15000, FWHM ~ 20 km s"^ at 
~ 1000 A) is more than one order of magnitude larger than that of the Chandra HRCS-LETG, 
and is sufficient to resolve the OVl2s-+2p absorption into at least two components (Figure 3a; 
Table 1). The total width of the OVl2s-.2p complex in the FUSE spectrum is FWHM=264±35 
km s~^, fully consistent with the upper limit on the width of the high ionization X-ray lines 
(Table 1). We also note that both position and width of the Ol3s_^3p(A = 1039.23) absorption 
line in the FUSE spectrum, arc clearly different from widths and positions of each of the two 
OVI components in the same spectrum (Fig. 3a,f), but again fully consistent with width 
and position of the 01(A = 23.489) line in the Chandra spectrum (Fig. 3f,c, Tabic 1). 

We then conclude that: (a) the high ionization lines, both in the UV and X-rays are 
likely to be produced by gas with similar dynamical properties, and broad range of dispersion 
velocities (from the UV profile), while (b) the low ionization atomic 01 lines (both in FUSE 
and Chandra) are produced by a slow-moving cloud of cold gas, with relatively low internal 
dispersion velocity, and so differs dramatically in both physical and dynamical properties 
from the high ionization gas. 

In the following we discuss only the highly ionized component. 
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Fig. 3. — Residuals, in velocity space, of the relevant portions of the FUSE (a,f) and HRCS-LETG 
(b,c,d,e,g) spectra of PKS 2155-304, to their best fitting continuum models. The residuals are centered, from 
the top to the bottom panels, on the rest-frame wavelengths of the resonant transitions: OVl2«^2p, OVIIxa, 
NelX^-a, OVIIIkc, OVII/f^, OI(A = 23.489). Bin widths are as follow: 2 km s^^ (FUSE: panel a,f), 190 km 
s~^ {Chandra: panels b, d, e), 160 km s~^ {Chandra: panel g), and 280 km s~^ {Chandra: panel c). 



3.2. Curve of Growth Analysis 



Absorption lines equivalent widths (EWs), and their ratios, can in principle provide 
invaluable information on the physical state of the gas, and allow one to distinguish between 
quite different scenarios. However these diagnostics can only be applied to unsaturated 
lines. Fortunately the UV OVl2s^2p hnes are clearly not saturated, since the ratio between 
the EWs of the doublet (Ai = 1031.9261 Aand A2 = 1037.6167 A), is fully consistent with the 
oscillator strengths (OSs) ratio of 2 (Savage et al, 2000). For the X-ray lines of OVII, OVIII 
and NelX, in principle, one could compare the EWs of the corresponding Ka and Kf3 lines. 
However, the errors on the EW s of the K(3 lines from these ions are too large. An alternative 
approach uses the Curve of Growth (CoG; e.g. Spitzer, 1978), computed for a given Doppler 
parameter, to see whether each line falls on the linear branch of the CoG. We ran our 
model for resonant absorption (Nicastro, Fiore and Matt, 1999) to produce CoGs for the 
OVl2s^2p(A = 1031.9261), the OVII^a, the OVIIIi^a doublet, and the NelX^a- We assume 
a Doppler parameter of b=200 km s^^, since this is the observed a of the OVl2s-->2p complex. 
The CoGs are plotted in Figure 4, along with the maximum possible measured EWs (best 
fit + 2(7) of each line (EWs are plotted in absolute value). For the OVl2s^2p(A = 1031.9261) 
complex we considered the total EW from the two best fit gaussians. All the lines EWs fall 
well within the linear branches of their CoGs. Hence the lines are not saturated, and the 
resulting column densities, listed in Table 1, should be reliable. ^ . 



4. Discussion 

The FUSE and Chandra spectra of PKS 2155-304 contain absorption features from local 
intervening ionized gas that span a broad ionization interval, from neutral (01), through 
mildly ionized (OVI, OVII), to highly ionized (OVII, OVIII and NelX). The neutral and 
the ionized components have very different dynamical and physical properties and so must 
belong to two different, spatially separated, systems. The profile of the OVI complex in 
the FUSE spectrum suggests that the ionized component is also made up of at least two 
different absorbers: one slow-moving, with relatively low internal dispersion velocity and 
possibly associated with a cloud in the Galactic disk, and one moving towards us with a 
velocity of about 150 km s""^ (~ 100 km s~^ transforming from the Heliocentric to the 



^The partial covering and scattering effects that can occur when the absorbing material is close to the 

source, as in AGN "warm absorbers" and broad absorption lines (e.g. Arav et al.. 2001), do not apply in 
this case, since PKS 2155-304 has z=0.116, while the absorbers are at their rest wavelengths, and so are 
physically well separated. 
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Fig. 4. — COGs for the detected UV and X-ray absorption lines: 0Yhs^2p{X = 1031.9261), OVII/fa, 
OVIIIko (doublet), and NelXft-Q,. A doppler parameter of b=200 km s~^ (the observed a of the OVl2s^2p 
complex) is assumed. Horizontal solid lines are the maximum possible observed EWs (best fit + 2a). 
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Galactic Standard of Rest frame ^ ) and with a similar internal dispersion velocity and 
possibly associated either with an ionized high velocity cloud in the Galactic halo, or with a 
WHIM filament. 

In the next sections we demonstrate that neither pure coUisional ionization nor pure 
photoionization models can account for the presence of strong OVII and NelX absorption 
with the observed intensity of the OVI lines, unless the three absorbers (one X-ray and the 
two UV) are associated with three different Galactic clouds with drastically different equi- 
librium conditions (with temperatures differing by an order of magnitude), and so locations. 
Alternatively we show that if there is one highly tenuous absorber, then the photoionization 
contribution by the diffuse X-ray background becomes important and can change the inter- 
nal ion abundance distribution to account for all the UV and X-ray observational constraints 
in a single absorbing medium. These solutions however, imply path lengths of the order of 
several Mpc, and so definitively locate the UV- X-ray absorber in the intergalactic medium 
(IGM), beyond our Galaxy on a scale of the Local Group, or somewhat bigger. 



For the photoionization contribution by the diffuse, extragalactic X-ray background, 
we adopt an EUV to 7-ray Spectral Energy Distribution [i.e. the specific fiux Fi,kg{E)] 
given by a broken power law with an exponential high energy cutoff (i.e. Fbkg{E) — 
[Xl(E)-°-/' +X2(E)-"''-'']e-(-^/-^/), with XI = for > E^, K2 = Q lox E < Eb, 
oisoft — 1-1) 0(hard = 0.4, Ef, = 0.7 kcV, and Ef = 50 keV; we also adopt a low energy cutoff 
at energies lower than the hydrogen ionization threshold Ehi)- We assume a normalization 
of 10 ph s"^ cm-2 keV"^ sr'^ at 1 keV (Parmar et al., 1999, Boldt, 1987, and Fabian & Bar- 
cons, 1992). Following the Hellsten et al. (1998) parameterization for the redshift evolution 
of the X-ray background, we can write the ratio between the diffuse background ionizing 
photons density and the electron density of a filaments at redshift z (i.e., the ionization 
parameter U) as: 



U = [{l+zf/inc\n;^ / diyin{Fbkg{iy)/hiy) ~ 1.4x10-^^(1+-^)^ = 0.7S-\nbh^/0.02)~\ 



where vhi is the frequency of the H ionization threshold. 

In eq. 1 we defined the local IGM overdensity compared with the mean density of the 



4.1. 



The XRB Contribution 




+ 0O 



(1) 



^http: / /nedwww.ipac.caltech.edu/forms/veLcorrection.htnil 
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Universe, Ue, SiS 5 — Ue /n;, where = (2 x 10-'^(1 + zyiQ^h'^/Om) cm-^ (Hellsten et al., 
1998). Hence, for typical ovcrdcnsitics of (5 = 10 (Hellsten et al., 1998, Dave et al., 2000), we 
expect, at any redshift, U 0.16^^, for fi^/i^ = 0.0125 (Walker et al., 1991; h = 0.5). Given 
the flat X-ray spectrum of the X-ray background, XRB photoionization greatly modifies the 
ionic distribution of a purely collisionally ionized gas at temperatures lower than logT ~ 6.4, 
for electron densities of rZe^lO"^ cm~^ (5;; 800). Photoionization is not important at typical 
ISM densities of 10~^ to 1 cm~^ (e.g. Spitzer, 1990; but see §4.5 for alternative ionizing 
sources) . 

4.2. Models 

We ran CLOUDY (vs. 90.04, Ferland, 1997) to build models for temperatures (in K) 
from logT = 5.1 — 8 and densities Ug = 10"'' — 1 cm"^. Figure 5 shows the ratios between 
the predicted ionic abundances of OVH/OVI (black curves), OVHI/OVI (green curves), 
OVHI/OVH (red curves) and NelX/OVII (blue curves). Sohd curves are for — 1 cm~^, 
while long-dashed curves correspond to — 10~^ cm~^ (i.e. 5 ~ 10). The thick sections of 
the lines delimit the 2a allowed intervals given the measured relative ion column densities 
derived from the corresponding EWs (see Table 1), using the formula - valid only when the 
lines are not saturated - EW{X^) oc A{X)nx^ (Nicastro, Fiore, & Matt, 1999), where A{X) 
is the relative abundance of the element X, compared to H, and nx^ is the relative density of 
the ion i of the element X. Figure 5 uses the total measured EW for the OVl2s^2p complex 
(i.e. narrow and broad components combined). 

According to the above formula, the ratio between equivalent widths of two lines from 
two different ions of the same element, then, depends only on the ratio of the corresponding 
equivalent widths. Instead, for ions from different elements, the relative ion densities are 
proportional to the inverse of their relative metallicity (i.e. nxi/nyi oc {EW{X^)/ EW{Y^))x 
{A[X) / A{Y))~^) . For the NelX/OVII ratio (blue lines) we show allowed regions for a solar 
Ne/0 ratio (magenta, dotted line, delimited by perpendicular tickmarks) and for 2.5 times 
solar ([Ne/0]5oior- = 0.158, Grevesse & Anders, 1989, Grevesse & Noels, 1993. Blue, dashed 
hne). This values are discussed in §4.3. 

4.3. Galactic versus Extragalactic Origin 

Let us consider first the Galactic "high-density" solutions, for which the gas is in pure 
coUisional ionization equilibrium. The ratio of OVII;<:q, to the entire OVl2s_»2p UV complex 
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LogT (K) 

Fig. 5. — Ratios between the predicted ionic abundances of OVII/OVI (black curves), OVIII/OVI (green 
curves), OVIII/OVII (red curves) and NelX/OVII (blue curves), for two electron densities: (a) Ue = 1 cm~^ 
(solid curves), and (b) rtg = 10^^ cm~'^ (long-dashed curves). Thick tracks on these curves delimit the 
2a allowed intervals for the measured relative abundances ratios for the same ions, derived for the total 
measured EW of the OVI2S— 2p complex. Ne/0 ratios of 1 (magenta, dotted line) and 2.5 (blue, dahed line) 
times solar are plotted. 
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(thick solid black curve) sets a quite stringent upper limit of logT < 5.7 (T in K) on the allowed 
temperature of the OVI-OVII absorber, if both the lines are produced by the same material. 
However, the OVIII and NelX ions are clearly too highly populated to be consistent with 
the intensities of the OVI and OVII lines, requiring logT> 6.1, and so no common solution 
exists for all the lines. We find that this is true for any ng^lO"^ cm~^. 

Let us now consider the extreme extragalactic "low-density" (ug = 10~^cm~^) solution 
(Fig. 5, thick long-dashed curves) where photoionization by the diffuse XRB overionizes the 
medium. These conditions allow the oxygen OVI, OVll, and OVlll ions to have a common 
temperature (shaded region in Fig. 5) for the relatively low temperature range 5.35 < logT 
< 5.50. A consistent solution for even the oxygen lines, requires the absorber to be highly 
tenuous (ne^5x 10~^ cm"^). A solar Ne/0 ratio (magenta, dotted intervals) is not consistent 
with these conditions, but a Ne/0 ratio of ~ 2.5 times solar (blue, dashed interval) provides 
consistency for all the ions (see §4.7). Increasing the Ne/0 metallicity ratio decreases the 
NelX/OVII relative density ratio estimated from the corresponding measured EW ratio - 
see §4.2 - and so lets the thick intervals in Figure 5 shift toward the left on the theoretical 
nNeix/novii curve, down to the region where a common solution for all the ions can be 
found. 

For a given temperature, we can derive the equivalent H column density of this gas. 
Assuming Ne/0 = 2.5 and logT = 5.45 gives Nh — 1.5 x 10^^[0 / H]^^^ cm^^. This translates, 
for a constant electron density of = 10~® cm~^, to a linear size of the absorbing cloud of 
gas of 4:.9[0 / H]q^ Mpc, consistent with a filament with an overdensity of ~ 10 located in 
the intergalactic medium near to our Galaxy. For O abundances of 0.3 x Solar (as found 
in clusters of galaxies, Sarazin, 1988) the size would increase to ~ 15 Mpc (~ 2 % of the 
distance to PKS 2155-304, and ~ 25 % of the distance to the first group of galaxies close to 
our line of sight to the blazar: Penton Stocke, & ShuU, 2000). 

4.4. Multi-Component Solutions 

The above discussion is based on the hypothesis that both the broad and narrow compo- 
nents of the OVI absorption complex are produced by a single, structured, absorber that is 
also producing the X-ray absorption. Here we explore alternative multi-component solutions. 

Let us suppose that only one of the two OVl2<i^2p absorption lines in the UV spectrum, 
is associated with the X-ray absorber. The broader OVI line contributes about 1/3 to the 
total observed equivalent width from the OVI system, while the remaining 2/3 are provided 
by the narrower fine (see Tab. 1). Based on the strengths of the OVl2s-+2p lines and on the 
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above discussion, it is then natural to associate the broader UV component with the X-ray 
absorber. In this case the narrow and more intense OVI component would be produced by 
a Galactic cloud of moderately ionized gas in collisional equilibrium. This "narrow" line gas 
must have a temperature of logT""'^'''°^J;5.0 if the size of the cloud is not to exceed 1 kpc, 
the size of a big halo cloud (e.g. Spitzer, 1990), given the measured EW and assuming a 
maximum Galactic ISM density of ~ 1 cm~^. The upper limit on the temperature, instead 
is set by the observed EW(OVI5,"::27/EW(OVIIxa) ratio, and this gives logT""''^'^ « 5.6. 

The second absorber should instead produce the broader and weaker OVI line, and all 
the X-ray lines. Wc tested this hypothesis by building a diagram similar to that in Fig. 5, but 
considering a variable fraction (from to 1) of the EW of the broad OVl2s^2p hue (Figure 6). 
We find that equilibrium, high density solutions exist only if a fraction of 5 to 25 % (at 2a) 
of the broad OVI fine is produced by the same gas producing the X-ray absorption. These 
solutions are limited to a range of temperatures around logT=6.4, whose width depends on 
the relative [Ne/0] abundance: A{logT) = 0.16 — 0.32 for [Ne/Ojo = 1 — 2 (dashed regions 
in Fig. 6). For logT = 6.4 and [Ne/0]© = 1, we find Nh = 2.2 x W^[0/H]-\ and so a 
linear size oi D — 7[0/H]q^ pc, for rze = 1 cm~^. 

Since only 25 % of the broad OVl2s-»2p can come from X-ray absorbing gas, a third 
ad hoc Galactic absorber similar in dynamical properties to the absorber imprinting the 
broad OVl2s^2p line, but with temperature differing by at least an order of magnitude from 
the temperature of the UV- X-ray absorber, is required to explain the remaining broad OVI 
absorption. Alternatively, a highly structured and, again, ad hoc multi-temperature plasma, 
with at least two zones separated by a sharp front, and with the same kinematics, would 
be needed to account for all the observed UV and X-ray lines. We can broaden the range 
of allowed temperature and increase the fraction of OVI absorption produced by the X-ray 
absorber, by allowing both the Ne/0 ratio to be supersolar, and the gas density to be lower 
than ~ 10~^ cm~^. However this converges to the extragalactic solution (§5.2). We find that 
an acceptable solution with a minimum size is obtained for logT = 5.7, ng = 6 x 10~® (i.e. 
5 ~ 60), and Ne/0 = 2, which gives Nh = Q x 1Q^^[0/H]o_1^q cm'^ and D = 3[0/H]q1^q 
Mpc. 

4.5. Alternative Photoionizing Sources 

As we have demonstrated (§4.1, 4.3) X-ray photoionization is needed to modify the 
ion relative abundance distribution from the pure collisionally ionized case (with logT = 
5.3-5.7), to allow for the simultaneous presence of OVI, OVII and OVIII, and also NelX. In 
order for the photoionization contribution to play a significant role, the ionization parameter 
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Fig. 6. — Fraction of broad OVI line (f(OVI)) produced by the same gas producing the X-ray absorption, 
as a function of gas temperature. The diagram shows the two regions [for [Nc/O]© = 1 (red), and [Ne/0]© = 
2, (blue)], in the f(OVI)-logT space, for which common UV-Xray Galactic solutions do exist: at most 25 % 
of the broad OVl2s->2p can come from a Galactic cloud of X-ray absorbing gas 
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has to be of the order of C/ ~ 0.1 (§4.1). Since the ionization parameter. U, is simply the 
ratio between the number density of ionizing photons and the electron density of the gas, 
U increases as either the electron density in the gas decreases (as in the IGM case) or the 
distance between the source of ionizing photons and the gas cloud decreases. So high values 
of U might also be reached in gas with typical ISM densities (i.e. rig ~ 10~^ — 1 cm~^) 
illuminated by other, discrete and locally much brighter hard X-ray sources than the XRB. 
In this section we explore this second possibility. 



4-5.1. Isotropic Nearby Source 

Based on the needed value of U ~ 0.1 we can put constraints on the brightness of 
the ionizing source as seen by us, and on the angular separation between the source and 
the gas cloud. Let us consider an intervening cloud of gas located at a distance r from 
us, along a given line of sight, and an ionizing source, located at a distance D from us. 
Let d be the distance between the ionizing source and the gas cloud. We can write: (a) 
L — AnD'^Fobs', (b) L — AncPFion', and (c) rtgO — d. In the above equations, L is the 
source luminosity, and Fobs and Fjo„ are the source fluxes as seen by us and the gas cloud 
respectively. These two fluxes are related to each other through the source-gas angular 
separation 9: Fion = {sin9)~'^Fobs- The ionization parameter U at the illuminated face of the 
cloud can then be written as a function of the source flux as seen by us: U = qn~^ {sinO)^'^ , 
where q = dEE^^ fobs{E) is the total flux of ionizing photons {Fobs = / dEfobs{E)). To 
maximize the contribution from hard X-ray ionizing photons we use as the ionizing continuum 
a cutoff power law with F = 1.5 from the H ionization threshold up to 100 keV, and F = 4 
from 100 keV up to infinity. This gives a total fiux of ionizing photons of ~ 1.3 x lO^F^^^"^*^ 
ph s~^ (where F'^^^^^ is the observed fiux between 2 and 10 keV, and we assumed Fgbs — 
3 X Flbs^^)- Our value U ~ 0.1 then gives: ?7 ~ 4 x l{)-'^n-^{Fl~^^^)cGs{sin-'^e) ~ 0.1. 
For ^ = 30 arcmin and gas densities in the ISM range = 10~^ — 1 cm~^, then, the 2-10 
keV observed flux must be as large as F^^^^ ~ 2 x 10~^ — 10"'' erg s~^ cm~^ (~ 10^ — 10^ 
Crab), to cflicicntly ionize the gas. For a 30 times smaller separation, this fluxes reduce to 
2 X 10^^ — 10^^ erg s^^ cm^^ (~ 0.1 — 10 Crab). No such source is currently seen. 

For halo distances of 30 kpc the implied 2-10 keV luminosities are L'^^^^iO = 30') = 2 x 
(10"^^ — lO'^'^) erg s~^ {d = 260 pc), far too large for any known stellar-size Galactic source, 
and L1^^^{9 = 1') = 2 x (10^^ — 10^°) erg s~^ {d = 8 pc), just about the maximum observed 
luminosity for the most luminous known Ultra-Luminous- X-ray sources (ULXs) or beamed 
microquasars. For the much closer distances of Galactic disk objects, say, 10 pc, these 
luminosities reduce to more normal values for X-ray binaries (i.e. L'^obs^i^ ~ ^0') = 2 x 
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(10^4 - 10=^6) erg s-\ and = 1') = 2 x (10^^ - 10=^^) erg s-^) but the physical 

separations between the source and the gas cloud become very small: d — 0.09 pc and 
d = 0.003 pc respectively. So, either very luminous sources, or very close source-to-gas 
systems are needed to efficently ionize a tenuous gas cloud in our Galaxy to the observed 
degree. 

We searched the ROSAT-WGACAT (White, Giommi, & Angelini, 1994) for bright X- 
ray sources in a 30 arcmin radius cylinder, centered on the line of sight to PKS 2155-304. 
We found 33 unidentified sources plus an FO star with ROSAT-PSPC count rates between 
0.0019 and 0.024 ct s~^, and angular distances from the line of sight from = 2.1 to = 20.4 
arcmin. The closest of these sources to the PKS 2155-304 line of sight has a count rate of 
0.0073 (source A: IWGA J2158.8-3011), while the brightest (0.025 ct s"^) is at a distance 
of 15.1 arcmin (source B: IWGA J2158.2-3000). We used PIMMS ^ to convert these count 
rates, into flux. To maximize the possible contribution at energies E> 1 keV, we assumed a 
hard X-ray continuum (a powerlaw with F = 1.5) heavily absorbed by a column of Nh = 10^^ 
cm~^ intrinsic neutral gas. This gave 2-10 keV fluxes of 10~^^ erg s~^ cm~^ (5 x 10~^ Crab), 
for source A, and 3.6 x 10~^^ erg s~^ cm^^ (1.8 x 10""^ Crab) for source B, at least (i.e. using 
rig = 10~^ cm~^) 4 and 5 orders of magnitude, respectively, lower than the values required 
to photoionize the gas at these angular separations. 

4-5.2. Beamed Nearby Source 

Let us consider a microquasar highly beamed into the direction of the gas cloud, and 
whose line of sight to us is somewhat perpendicular to the beam direction. The probability 
for a microquasar to be beamed in the direction of the gas cloud, and located within a 20-130 
pc radius sphere from the gas cloud, depends on the distribution of the parent poulation, 
that, for microquasars, is largely unknown. We can compute here the minimum beaming 
factor needed, in the above optimal geometrical configuration, for source A or B to efficiently 
ionize a cloud of gas along the hue os sight to PKS 2155-304. The luminosity amplification 
due to beaming, in the direction of the jet, is given by: Lbeamed — Lintr x i'^lT- The 
observed luminosity along a perpendicular direction is insted de-boosted by the factor: Lobs — 
Lintr X (Urry & Shafer, 1984). In the above equations, Lintr is the source intrinsic 
luminosity, 7 is the Lorentzian factor, and p ~ 4 (Urry & Shafer, 1984). Solving for the 
intrinsic luminosity, and using as maximum value for the beamed luminosity of a microquasar 
Lhcam = 5 X lO^*' erg s^^ and the observed value of Lobs = 4 x 10'^^ erg s^-^ (for source B, at 



®http://legacy.gsfc. nasa.gov/Tools/w3pimms. html 
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30 kpc), we obtain a minimum relativistic factor of 7 > 3, to compare with the estimated 
upper hmit of 7 = 3 for the microquasar GRS 1915+105 (Mirabel & Rodriguez, 1999). Any 
deviation from the above ideal geometry, would give a much larger Lorentz factor. 

4.5.3. Historical Nearby Source 

Low density gas has long recombination times. A source may thus have been much more 
luminous in the past and ionized the gas then, from which the gas has not yet recovered. Let 

us first suppose that the most luminous X-ray source we are currently seeing in the ROSAT 
archive is a binary system that went through a ULX phase (King et al., 2001) in the past 
(from the calculation in §4.5.1, for source A and B we need L(outburst)/L(now)> 10^ - 10^ 
respectively ^ ). The recombination time for the OVll-OVlll ions, in gas with electron 
density of rig = 1 cm~^ is of about t 

rec ~ 10^ years (Nicastro et al., 1999). Typical ULX 
phases are thought to last for ~ 10^ years, and in our Galaxy the total estimated rate of 
ULXs is of 1 every 10^ years (King et al., 2001). So, allowing for a maximum delay of 
tdeiay — Uec — 10^ ycars, wc may expect N = 2 ULXs in all the Galaxy over a period of 
— tphase + tdelay = 2 X 10^ ycars. Assumiug, for the halo, a 30 kpc radius sphere, this 
gives a ULX density, over the same period, of n = 2/((4/3)7r x (30)%pc) ~ 1.8 x 10~^ kpc~^. 
So, the maximum number of ULXs expected over a period At in a random Rmax = 20 — 130 
pc radius sphere in our galaxy (the physical distances between sources A and B and the line 
of sight to PKS 2155-304, if both sources are put at 30 kpc from us, in the outskirts of the 
Galactic halo) is of Nmax = 6 x 10~^° — 2 x 10^'', making prior ionization by a powerful 
transient source highly unlikely, albeit with a posteriori statistics. 

4.5.4. Historical AGN at Galactic Center 

A final possibility is that in the past the center of our Galaxy went through an AGN 
phase, producing far higher luminosities and photoionizing any low-density {n^ ~ 10~^ cm~^) 
clouds in the Galactic halo. At these densities we have trec{OVII, OVIII) ~ 10^ years. Let 
us then suppose that less than 10^ years ago the 2 x 10^ Mq black hole in the center of our 
Galaxy (Chakrabarty & Saha, 2001) was emitting at its Eddington luminosity of 2 x 10^^ 
erg s"-*^. Assuming a typical Spectral Energy Distribution (SED) for AGN (Mathews & 
Ferland, 1987) this would produce an ionization parameter of Uagn — 2.4 x 10^^ x at 



''Such outbursts have been observed in X-ray novae between quiescent and active phases, e.g. Narayan, 
Garcia, & McClintock, 2001. 
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the surface of Ug — 10~^ cm~^ halo clouds located at a distance D from the central AGN. 
To produce the observed line ratios, the photoionization contribution must be such that 
Uagn ~ 0.1 at the illuminated face of such a cloud of gas. This gives a distance for the gas 
cloud of D ~ 36 kpc, consistent with the outermost size of the Galactic halo. This solution 
requires that all halo clouds would show UV and X-ray absorption with similar physical 
properties as those observed in PKS 2155-304. Using the much lower luminosities estimated 
by ASGA for the central AGN in our Galaxy about lO'' years ago (i.e. L ~ lO''^ - 10^^ 
erg s~^, Koyama et al., 1996), gives distances of 1-5 kpc, well within the Galaxy's halo, 
and therefore not consistent with the the galactic coordinates of PKS 2155-304 (given our 
distance of about 10 kpc from the Galactic center). 



4.6. Properties of the WHIM Gas 

Prom our analysis of the combined UV and X-ray data, we conclude that the X-ray 
absorption and at least the broad OVl2s_»2p are likely to be produced in tenuous extragalactic 
medium, with rig ~ 6 x 10~^ cm~^. Prom the position of the OVl^f_^2p assuming that 
the feature lies within a few Mpc from our Galaxy, the gas is falling toward our Galaxy with 
a bulk velocity of ~ 100 km s~^ (in the Galactic Standard of Rest frame). The feature has a 
dispersion velocity (~ 300 km s~^) comparable with the bulk velocity, and a temperature of 
logT~ 5.7. The implied overdensity is 5 ~ 60 and the total linear size of the feature along 
our line of sight is about '^SqqIO / H]q Iq Mpc. 

The Ne/0 abundance can be reduced to Solar if we allow for inhomogeneities in the flow, 
since inhomogeneities increase the upper limits on both the temperature and the density, 
and so also decrease the lower limit on the linear size of the feature (from 3 [O / H]q1^q Mpc 
to few X lO^[O/if]o3x0 kpc). Such inhomogeneities may arise near our Galaxy halo, where 
the gas may mix with denser medium, cooling and producing part of the observed OVl25^2p 
absorption. 

The above overdensities and linear sizes are consistent with those predicted for the 
WHIM by hydrodynamical simulations for the formation of structure in the local Universe 
(e.g. Hellsten, 1998), and in particular in our Local Group (Kravtsov, Klypin, and Hoffman, 
2001). According to these simulations, the local IGM is concentrated in filaments with 
overdensities of (5 ~ 5 — 100 connecting already virialized structured (galaxies and clusters 
of galaxies), that were shock- heated during the collapse of the density perturbations. The 
average sizes of these perturbations are of the order of 30-50 Mpc in the local Universe, but 
the total extent of the density peaks are usually not larger than a few Mpc. The velocities of 
the matter in these filaments, should be of the order of the proper motions of the virialized 
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structures that they connect, ~ 100 — 200 km s ^. We note that these velocities are of the 
same order as those observed on the blue side of the structured and complex OVl2s_»2p profile 
in the UV spectrum of PKS 2155-304. 

We also note that both the dynamics of this gas (as derived from the OVI line) and 
its physical state are very similar to those found (using UV absorption hues only) for the 

intervening ionized absorber along the line of sight to the quasar H 1821+643 (Tripp ct al., 
2001): internal dispersion velocities of the order of ~ 150 km s~^ (FWHM of the broad Lya 
line), temperature in the range 5.3 < logT < 5.6 (assuming pure collisional ionization), and 
pathlcngth of about 3 Mpc (using Ho = 50 km s~^ Mpc~^). We predict that the planned 
high resolution and high quality X-ray spectra of this object, should then show OVlli^-Q 
absorption with an intensity similar to that found here. 



4.7. The Ne/O problem: Dust Depletion 

In virtually all the cases discussed, the Ne/O ratio in the absorber along the line of 
sight towards PKS 2155-304 is required to be higher than solar. There arc various ways to 
reconcile with this problem: (1) Type II supernovae give rise to Ne/O ratios much higher 
than Type I supernovae, compared to Si (Tsuro et al., 1997 and references there in); (2) O can 
be depleted on to dust grains, but noble gas like Ne would not. If superwinds from galaxies 
are responsible for the enrichment of the IGM, then in order to produce high gaseous Ne/O 
ratios compared to solar, they must be able to produce a considerable amount of dust before 
leaving the high density galaxy environment. (3) In the modeling above, we have assumed a 
simple parameterization of ionizing background. The true shape of the background is clearly 
more complex, and a harder background would result in a large Ne/O ratio. Any or all of 
these factors might be playing a role in the observed value of [Ne/O]. 

It is important to notice that dust can easily survive at the kinetic temperatures (i.e. 
few X 10^ K) and densities (i.e. ~ few x 10^^ cm^^) of such a warm IGM filament. In such 
regime, the sputtering lifetime is inversely proportional to the gas density, and largely inde- 
pendent from the dust grain composition and gas temperature (Draine & Salpeter, 1979a,b). 
The survival time of a small interstellar grain of size ~ 0.01 /im will be larger than 10^° yr 
in the IGM filament, and even larger for grains of larger size. We estimated the reddening 
expected along the line of sight to PKS 2155-304, based on (a) the observed equivalent H 
column density, (b) the typical size of a silicate dust grain, and (c) a typical ISM gas/dust 
ratio. In this condition the estimated E(B-V) is of about 2 x 10~^ mag, fully consistent with 
the observations (Schlegel et al., 1998). 
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Finally we note that [Ne/0] overabundances of ~ 2 — 3, compared to solar, have already 
been observed in several astrophysical environments (e.g. Paerels et al., 2001, Brinkman 
et al., 2001), strongly suggesting that the non-O-depleted values of this quantity, measured 
in stelar atmospheres, are not representative of the average gaseous [Ne/0] ratio in the 
Universe. 



5. Conclusions 

In this paper wc report the first discovery of the "X-ray Forest" of absorption lines pro- 
duced by highly ionized intcrgalactic medium. High significance X-ray absorption lines of 
OVII, OVIII and NelX are detected along the line of sight to the bright blazar PKS 2155-304, 
and are associated with a known OVl UV absorber. Wc demonstrate that the dynamical 
properties of the X-ray and the UV absorber are fully consistent with each other, and that a 
reasonable common and self-consistent physical solution can be found only if photoionization 
and coUisional ionization both contribute to the ionization of the absorbing gas. This requires 
electron densities of about 6 x 10~^ cm~^ for the diffuse X-ray background to be significant. 
This low denisty requires a linear size, along the line of sight, of the order of ^[0/H]q1^q 
Mpc. This clearly locates the absorber outside our Galaxy in intcrgalactic space. We demon- 
strate that both the dynamical and physical properties of such an absorber are remarkably 
consistent with those predicted for the low redshift warm phase of the IGM as predicted by 
hydrodynamical simulations for the formation of structures in the Universe. Finally we find 
that solutions with a Ne/0 ratio of about 2.5 times solar are favored, suggesting type-II 
Supernova enrichment, or the presence of dust in the IGM. 
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